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Abstract 

Valence energies for crystalline C, Si, Ge, and Sn with diamond structure 
have been determined using an ab-initio approach based on information from 
cluster calculations. Correlation contributions, in particular, have been eval- 
uated in the coupled electron pair approximation (CEPA), by means of incre- 
ments obtained for localized bond orbitals and for pairs and triples of such 
bonds. Combining these results with corresponding Hartree-Fock (HF) data, 
we recover about 95 % of the experimental cohesive energies. Lattice con- 
stants are overestimated at the HF level by about 1.5 %; correlation effects 
reduce these deviations to values which are within the error bounds of this 
method. A similar behavior is found for the bulk modulus: the HF values 
which are significantly too high are reduced by correlation effects to ~ 97% 

of the experimental values. 
71.45.Nt, 31.20.Tz, 71.55. Cn 
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I. INTRODUCTION 



Several methods are presently available for performing ab initio calculations for solids. 
Most frequently used is the density-functional method, with a local-density approximation 
(LDA) for the exchange and correlation contributions to the total energy§. This method 
yields good results for solid-state properties, on average, but a systematic improvement 
is difficult. In the past few years it has become possible to determine self-consistent-field 
(SCF) calculations for solids using the exact non-local Hartree-Fock (HF) exchangei. While 
the results are often inferior to LDA because electron correlations are neglected, it has 
the advantage of yielding a well-defined mean-field wave-function, which can be used as 
a starting point for treating many-body effects. In the Quantum Monte Carlo (QMC) 
approach!, one multiplies a mean field function (Slater determinante) by a Jastrow factor 
which explicitly introduces inter-electronic coordinates. Another approach of incorporating 



wave-function are used to admix suitable one- and two-particle excitations to the mean- 
field ground state. A method closely related to the ideas of the LA is the incremental 
expansion of the correlation energy; here information on local excitations in clusters is made 
use of, which are accessible to an accurate quantum-chemical configuration interaction (CI) 
treatment. This way ~ 85 % of the correlation contribution to the cohesive energy E co h has 
been recovered for diamond! and crystalline silicon!. In the present paper we extend the 
application of the method of increments to two more systems - germanium and grey tin - 
and to other properties beside E co h. we determine, for all of the group IV semiconductors, 
the influence of electron correlation on the lattice constant a and the bulk modulus B. 

In order to discuss the influence of correlation, we need reliable HF values for the ground 
state properties of these solids. We supplement the results of the Torino group0 by cluster 
calculations which are described in Section 2. In Section 3, we sketch the method of in- 
crements and report on computational details of the correlation treatment. The results for 
group IV semiconductors are discussed and compared with experiment and other calculations 



electron correlation is the local ansatz 




where local operators acting on the SCF 
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in Section 4. Conclusions follow in Section 5. 



II. HARTREE-FOCK CALCULATIONS 

Obtaining accurate HF energies for solids is still a difficult computational task. The 
Torino group of Pisani and co-workers developed a method which enables one to perform 
SCF-LCAO (linear combination of atomic orbitals) calculations for periodic systemJl. How- 



encountered, if one uses Gaussian functions with small exponents as normally contained 
in molecular basis sets; moreover, higher polarization functions than d are presently not 
available. 

Another possibility for obtaining information on solids are cluster calculations. They 
do not account for the infinite extension of the solid, but one can use standard quantum 
chemical program packages for achieving high accuracy, without any restriction with regard 
to the basis set. In this paper, we use a cluster method in conjunction with an energy 
partitioning approach which was developed for diamond^ and is here extended to silicon, 
germanium and grey tin. The systems considered in this scheme are fragments of the di- 
amond lattice, with dangling bonds saturated by hydrogen atoms. More specifically, one 
chooses clusters with closed structure where each of the X atoms (X=C, Si, Ge, Sn) has at 
least two neighbouring X atoms. The smallest one which meets this condition is X 6 H 12 , the 
biggest one which was accessible in our calculations is X35H36. (All clusters of the present 
SCF study are shown in Fig. 1.) The X-X distances were taken from the solid (rQQ=1.544 



(pvdz) basis set (9s4pld)/[3s2pld] of DunningpJ was selected. For Si, Ge and Sn we em- 
ployed 4-valence-electron pseudopotentials simulating the atomic cores, together with the 
corresponding optimized basis sets of double-zeta quality (dz) (4s4p)/[3s3p] for the valence 



ever, as a 



certain drawback of their program package 




convergency problems are 




3 



electronsEl; we added one d-function in each case whose exponent was optimized in CI calcu- 
lations for XH 4 and X 2 H 6 (Si: 0.40; Ge: 0.32; Sn: 0.23). For hydrogen we chose Dunning'sEO] 
double zeta basis (4s) /[2s], without polarization function. We calculated the total HF energy 
^total f° r eacn cluster, using the direct SCF program package Turbomole0 of Ahlrichs 
and co-workers. In order to obtain an estimate of that part of E^ ^ which can be at- 
tributed to a 'solid-like' X atom of the cluster, surrounded only by X atoms, we employed 
the energy partitioning method of Refs.0&0: -^total is approximated by a sum of energy 
contributions of the bare X atom, XH and XH 2 groups 

^total = n X E X + n XH^XH + n XR 2 E XR 2 W 

where (^XH' n XH 2 ) are ^ ne numDers of atoms in the cluster with zero (one, two) 
neighbouring H atoms. We determined the quantities Ey^, -ExtF and ExH 2 adjustment 
to the SCF results -Etotal °^ ^ ne three largest clusters (X22H28, X26H30 and X35H 36 ). We 
checked the resulting values for the group energies using the total SCF energies of the other 
clusters: the quantity cr in Table 1 is the (average) difference of these energies to the values 
evaluated from the group contributions according to (1). Finally, the HF cohesive energy of 
the solid per unit cell was calculated from 

^coh = n uc(^X _ E atom) ( 2 ) 

where -Eatom * s the SCF energy of the free atom determined at the same level as EV (i.e. 
using the same basis set), and nuc is the number of atoms per unit cell (2 for the diamond 
lattice). The results for E^^ of the group IV semiconductors are listed in Table 1. The 
error a gives an estimate for the finite size effect; it is 3 % for Sn and smaller for the other 
compounds. Hartree-Fock results from literature^, calculated with Crystal, are in good 
agreement (to ~ 1%) with our values for diamond and silicon (although smaller primitive 
basis sets were used in Ref. 0) , but there is a difference of 7 % for Ge. 



III. CORRELATION-ENERGY INCREMENTS 
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A. Formalism 



In the following we give a derivation of the incremental expansion for the correlation 
energy which supplements the one given in RefsJH; it is similar to the one proposed in 
Ref.El. The present derivation is more formal, but it shows clearly that the method applies 
to infinite solids and should not be considered as being merely related to calculations for 
finite systems. We use thereby the CEPA-0 (see e.g. Ref.S) which is particularly suitable 
for our purpose. 

We start from a Hamiltonian H which can be divided into two parts 

H = H + H 1 . (3) 

The ground state of H is supposed to be known; in our case it is the HF ground state with 
the corresponding wave-function $ - We define a product of two operators A and B in the 
Liouville space as follows: 

(A\B) = (%\A^B\%) C = (A^B) C (4) 

The upper script c indicates that the cumulant of the expectation value is taken, which is 
given by 

(A) c = (A) (5) 
(AB) C = (AB) - (A)(B) (6) 
etc. 

For further details see RefM By using we may write the exact ground-state energy E 
in the following way: 

E=(H\n) = E RF + (H 1 \n), (7) 

where Q plays the role of the wave operator which describes the transformation from the 
HF ground state to the exact ground state. 
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If we have a solid with well-defined bonds, we can express the HF ground state $o i n 
terms of localized orbitals and label those orbitals by a bond index i. We define operators 
Ai, where i should being consider as a compact index which includes the bond i as well as 
the one and two particle excitations of bond i, and Aij, which describes the two particle 
excitations where one excitation is out of bond % while the other is out of bond j. Within 
the restricted operator subspace spanned by Ai and A^ we make for Q an ansatz of the form 

\Q) = \1 + J2 n i A i +J2 n iJ A ij)- ( 8 ) 

i ij 

This choice represents the coupled electron pair approximation at level zero (CEPA-O)il. 
The parameters rij and determined from the set of equations (Ai\HQ) = and 

(Aij\HQ) = 0@. With (D this implies 

= (A k \H) + Y,n l (A k \HA l ) + Y,n lj (A k \HA ij ) 

i ij 

= {A kl \H) + Y,n l {A kl \HA l ) + Y,n l3 {Aki\HA l3 ). (9) 



The method of increments provides a scheme, in which this set of equations and hence the 
correlation energy is evaluated in a hierarchical order. 

a) First all electrons are kept frozen except for the ones e.g. in bond i. The operators Ai 
describe the corresponding excitations of these two electrons and Eq. @ reduce to 

0=(^|fT)+n i (1) (A-|^)- (10) 

Within this approximation the are independent of each other and the correlation 
energy becomes 

E COYY = ^2 e i (11) 

i 

with 

ei = nf\H x \Ai). (12) 
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b) In the next step we correlate the electrons in two bonds, e.g. i and j. The correspond- 
ing rS r) parameters are determined from the coupled equations 

= (A\H) + n^iAlHAi) + nfiMHAj) + nf^A^HA^) 

= (A\H) + nf\A 3 \HA t ) + nf (A^HAA + ^{A \HA tJ ) 

= (A tJ \H) + n^iAjlHAi) + nfiA^HAj) + nf^A^HA^) (13) 

Again, the increments 5rii = rip — rif 1 and 5rij = nf — nf 1 are treated as independent 
of each other in this approximation, and we have 

4o } rr = E* + ^£^y (14) 
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where 



Ae ij = e ij -(e i + e j ). (15) 



and 



Si 



H^A + nfAj + n^A^). (16) 



c) Analogously we calculate the three bond energy increment, which is defined as 

Ae ljk = e ljk - (Aq + Ae j + Ae k ) - (A Qj + Ae jk + Ae lk ). (17) 

The correlation energy €ij k is that obtained when all electrons are kept frozen except 
those in bond i,j and k. Again, the increments Ae^ are treated as being independent 
of each other. 

The total correlation energy within this approximation is the sum of all increments. 

= E e i + ^E A ^ + 7^ E A^- fc + ... (18) 

i ij ijk 
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It is obvious that by calculating higher and higher increments the exact correlation energy 
within CEPA-0 is obtained. 

The method just described is only useful if the incremental expansion is well convergent, 
i.e. if increments up to, say, triples are sufficient, and if increments become rapidly small with 
increasing distance between localized orbitals. These conditions were shown to be well met 
in the case of diamond and siliconlli, but have to be checked again for germanium and grey 
tin here. Ideally the increments should be local entities not sensitive to the surroundings. 
We use this property to calculate the energy increments in finite clusters. If they can be 
proven to be well transferable even between finite clusters, such cluster calculations may be 
extrapolated to the corresponding solid state quantities. 



B. Computational Details 

In this section we give computational details characterizing our ab initio calculations for 
the correlation energies of group IV semiconductors with diamond structure. 

First we select suitable fragments of the diamond lattice, and, as in the SCF calculations, 
we saturate the dangling bonds with hydrogen. We have to take much smaller clusters, how- 
ever, for the correlation treatment (see Fig. 2) than in the SCF case, because CI calculations 
for clusters as large as X 35 H 36 would be prohibitive. 

As a second step we perform, for each cluster, a standard SCF calculation (using the 
program package Molpro@) and localize the bond orbitals, applying the Foster-Boys 
criterion^, within the occupied valence space in Ci symmetry. Following that we calcu- 
late CEPA-0 energy increments by successively correlating more and more of the localized 
X-X bond orbitals (LMO) as described in the previous subsection. 

Two different basis sets are used: Basis set A is the same as in the SCF calculations. 
An extended basis set B has been generated by replacing the single d function of basis set 
A by a 2dlf polarization set with an optimized / exponent (C: d: 1.0970, 0.318; f: 0.76; Si: 
d: 0.8, 0.23; f: 0.35; Ge: d: 0.6, 0.15; f: 0.43; Sn: d: 0.45, 0.12; f: 0.30). In order to check 



8 



the quality of the basis sets we performed test calculations for the one-bond increment Ae.j 
of the Ge 2 H 6 cluster (Table 2). One sees that the sp basis set is well chosen; enlarging the d 
space and adding an additional / function yields an enhancement of about 25% , supplying 
a polarisation set 3d2flg@ adds another 20% . 

Next we have to check the transferability and the convergence of the increments. We will 
discuss them for Sn since this is the most critical case. The correlation energy increments 
obtained for Sn are listed in Table 3 together with the weight factors appropriate for the 
solid. We observe the following: 

a) The convergence of the incremental expansion is quite rapid both with respect to the 
number of the bonds involved (e.g. cluster 4: e% ~ 2.5Ae23 ~ 300Aei23 ~ 1500Aei234) 
and with respect to the distance between the bonds (e.g. cluster 4: Aei2 ~ 10Aei3 w 
50Aei4). In view of these findings, we restrict the energy increments up to the third 
order for adjacent bonds and up to the second order for next nearest through third 
nearest neighbours. 

b) The transferability of the increments is reasonably good. There are changes of 1 • 10~ 3 
a.u. for one-bond increments between the smallest (Sn 2 H 6 ) and the largest (Sn 6 H 14 ) 
clusters considered. For the most important two-bond increment the difference be- 
tween Sn 3 H 8 and Sn 6 Hi 4 is about 1 ■ 1CT 4 a.u.. As these changes are of different sign, 
the total effect on the cohesive energy of the solid is not larger than 3 • 10~ 3 a.u. per 
unit cell. 

Finally, in order to test the quality of the CEPA-0 results, we performed calculations for 
GesHg at different levels of correlation treatment (Table 4). We tried CEPA-1 and CEPA-2 
as well as the coupled cluster method with single and double excitations (CCSD) and even 
with additional triple excitations included in a perturbative way (CCSD(T)). The effect on 
the one-bond increment is of the order of 3 • 10~ 4 a.u., and ~ 2 ■ 10~ 4 a.u. for the two-bond 
increment. Thus, the errors are of the same order of magnitude as those caused by lack of 
transferability. 
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It is clear, from these remarks, that the accuracy of our method has its limitations, in 
practical applications. On the other hand, there are formal arguments why errors due to the 
truncation of the incremental expansion should be small: higher than two-bond increments 
involve triple excitations which do not directly couple to the Hartree-Fock ground state; 
two-bond increments between localized orbitals at large distances are of the van-der-Waals 
type, with a rapid decrease as \. Moreover, there are means for numerically controlling 
these errors: the \ law just mentioned, e.g., lends itself to an approximate summation of 
neglected two-bond increments; energy variations between different finite clusters allow for 
an easy estimate of transferability defects of individual local increments. 

Summarizing these considerations, one may conclude that the greatest remaining error 
of the results, to be discussed in the next section, is due to limitation of the one-particle 
basis sets. 



IV. RESULTS AND DISCUSSION 

Applying the method of increments as described in the preceding section, we have de- 
termined correlation contributions to cohesive energies, lattice constants, and bulk moduli 
for all of the group IV semiconductors. The increments were always taken from the largest 
possible cluster (cf. Table 3 and Fig. 2). 

The correlation contributions to the cohesive energies were obtained as E^°j^ = E^®^ ~ 
2E^^ n per unit cell of the diamond lattice. The results for the two different basis sets are 
shown in Table 5. For basis set A we obtain about 65 % of the 'experimental' correlation 
energies (defined here as the differences between the experimental cohesive energies and the 
corresponding HF values of Sect. 2). The larger basis set B yields a substantial improvement, 
to about 85 % . Combining the HF results with the correlation contributions we recover 
about 95 % of the experimental cohesive energies^. (The experimental values in Table 5 has 
been corrected for the phonon zero point energies f^g©!) derived from the Debye modelill 
(C: 1860 K; Si: 625 K; Ge: 360 K; Sn: 260 K)fi) 
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For comparison, we have also listed in Table 5 results from the literature which have been 
obtained with other methods. The Local Ansatz (LA) which also uses a CEPA-0 scheme 
yields smaller cohesive energies for all compounds.^ (The LA values in Table 5 have been 
obtained by adding the correlation contributions to E co ^ calculated in ReH to the SCF 
values of the present work (Sect. 2).) Our basis set A is comparable with the basis used in 
LA, but still the correlation contributions differ by about 1 eV for Si, Ge and Sn. lda!@ 
overestimates the cohesive energies by ~15 %. The QMC result! for diamond is excellent 
and also very good for silicon (with an accuracy of 4%); note, however, that in both cases 
the HF cohesive energy of Ref.l is lower by ~ 1 eV than that of Ref.t3 and of this work. 

In the next step we evaluate lattice constants (Table 6). They have been determined by 
varying the X-X distances both in the HF and in the CEPA-0 calculations. More specifically, 
to obtain the HF lattice constant we varied all X-X distances of the largest cluster X35H36 
and minimized the SCF energy with respect to the interatomic distance. At the HF level, the 
lattice constant is larger than in the experiment, except for diamond. Our results are in fair 
agreement here with the Crystal SCF calculations by Causa and ZupanS: the deviations 
are ~ 0.5% for C and Si, but ~ 1% for Ge, see Table 6. By using pseudopotentials for the core 
electrons in the present work, any effects of the latter beyond the frozen-core approximation 
are excluded from the outset. However, calculations for small molecules show that core 
polarization is very important for bond lengths of Ge and Sn compounds^. We simulated, 
therefore, this effect using a core polarization potential (CPP)§ and studied the influence on 
the X-X distance for different clusters. We find that the changes are nearly independent of 
the cluster size. By transfering these changes to the solid we obtain the results listed in row 
b) of Table 6. The influence of core polarization is seen to significantly increase within the 
group Si, Ge and Sn. Valence correlations affect the lattice constant in two different ways. 
Correlation contributions calculated with a minimal basis set ('inter- atomic correlations') 
enlarge the lattice constant. The increase for diamond is mostly due to them. Intra-atomic 
correlations, on the other hand, decrease the lattice constant; extended basis sets with many 
polarization functions are needed to describe them accurately. The total influence of valence 
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correlation seems to become less important when going from Si to Ge and Sn. The final 
values are within the error bounds of the pseudopotential and CPP (~ 0.02 A), and within 
the error of the limited basis set which can be roughly estimated from the difference between 
basis sets A and B to be < 0.02 A. The LA overestimates the lattice constant, since with 
the small basis set used in Ref.il the intra-atomic correlations could not be described very 
well. The LDA yields lattice constants which are systematically too small by about 1%. The 
QMC result is again excellent for diamond, slightly less good for silicon, but still reaching 
the experimental value within their error bounds. 

As a last property, we consider the bulk modulus B = Vj^, which describes the response 
of the solid to a homogeneous pressure. For the diamond structure one can easily derive the 
following expression 

B = (i& ~ £ »hM- da) 

We evaluate the bulk modulus at the experimental lattice constant a, so that the second 
term in flT9| ) is small but not zero. We obtain the HF bulk modulus again from the largest 
possible cluster. Compared with experimental (see Table 7) all values are too high. The 
core-polarization effect has been taken from smaller clusters, which causes an uncertainty 
of about 3 %, but the total effect of the CPP is much greater so that one can justify this 
approximation: for Sn, it reduces B by 14 %. Valence correlation reduces the bulk modulus, 
too, especially for Ge. Overall, we obtain results for the bulk modulus which are a little too 
small but still within few percents of the experimental data. 

V. CONCLUSIONS 

We have determined ground state properties (cohesive energy, lattice constant and bulk 
modulus) of group IV semiconductors, both at the HF and the CEPA-0 levels. The HF 
results have been obtained using a cluster method and an energy-partitioning ansatz which 
works well for solids with covalent bonds as those considered here. Electronic correlations at 
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the CEPA-0 level are described with the method of local intra- and inter-bond increments, 
which allows for a systematic improvement of accuracy towards the fully correlated solid- 
state limit. The results show that this method works well for all solids with diamond 
structure, with the transferability only slightly deteriorating down the fourth column of the 
Periodic Table. Not only the cohesive energy E^^, but also the lattice constant a, and the 
bulk modulus B, have been calculated with quantum chemical accuracy, to about 5% for 
E co /i, 0.5% for a, and 3% for B. Work is underway in our laboratory to apply the methods 
presented in this paper to solids with zinkblende structure, such as GaAs. 
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TABLES 

TABLE I. Hartree-Fock cohesive energies per unit cell (in a.u.), with average deviations a in 
parentheses (cf. text) 





C 


Si 


Ge 


Sn 


this work 


-0.3947 


-0.2273 


-0.1560 


-0.1323 




(0.0036) 


(0.0032) 


(0.0034) 


(0.0044) 


Refill 


-0.3984 


-0.2253 


-0.1672 





TABLE II. Test calculations for the Ge basis set: the SCF energy, EgQp, and the one-bond 
CEPA-0 correlation-energy increment, Ae^, for Ge2Hg (in a.u.) 



basis set 


^SCF 




A 


(4s4p)/[3s3p] 


Id 


-10.831779 


-0.01845 


B 


(4s4p)/[3s3p] 


2dlf 


-10.836493 


-0.02265 




(4s4p) 


2dlf 


-10.836633 


-0.02301 




(6s6p) 


2dlf 


-10.838324 


-0.02317 




(6s6p) 


3d2flg 


-10.841788 


-0.02750 
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TABLE III. Correlation-energy increments for Sn (in a.u), determined at the CEPA-0 level 
using basis set A. For the numbering of the source clusters and bonds involved, see Fig. 2. 





Son tcp rlnstpr/ 


Increment 


Wipiffht, factor 




UU11LL UIUlLdlb 




1UI lilt! bUHU. 




1 /1 


-U.UloOZD 


/I 




9 /1 
2/1 


-U.UloZO 1 






A /9 

4/2 


-u.ui / yoo 






8/1 


n ni 7/i 7r 

-U.UI / 4 / 








-U.UUO ( »4 






4/2,3 


-U.UUb f 38 






8/ 1,2 


-u.uuoyuo 






4/1,3 


n nnn7ni 
-U.UUU lUL 


1 o 

Iz 




8/2,5 


-U.UUU 1 15 






7/1,3 


-U.UUU5 ID 






8/2,4 


-U.U0U553 






4/1 4 


n nnni49 


1 9 




£ /1 /I 
5/ 1,4 


n nnm ha 

-U.UUU1U4 


/I £ 




f; /1 /i 
0/1,4 


n nnm 1 1 

-U.UUUl 1 1 


9/1 




7/1,4 


-U.U0U194 


12 




8/1,2,3 


U.U00475 


8 




8/1,2,5 


-U.U0U038 


12 




4/1,2,3 


-U.U0U055 






8/1,2,4 


U.U0U055 


24 




4/1,2,3,4 


-U.U0U011 
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TABLE IV. Test calculations for the one-bond increment and the two-bond increment Ae^ 
between nearest neighbours, at various levels of correlation treatment, for GesHs (a.u.) 







Af 

^3 


CEPA-0 


-0.01824 


-0.00651 


CEPA-1 


-0.01798 


-0.00631 


CEPA-2 


-0.01798 


-0.00650 


CCSD 


-0.01798 


-0.00613 


CCSD(T) 


-0.01798 


-0.00645 



TABLE V. Cohesive energies per unit cell (in a.u.); deviations from experimental values (in 
percent) are given in parentheses 





C 


Si 


Ge 


Sn 


this work, basis set A 


-0.5077 


-0.2996 


-0.2346 


-0.1979 




(92 % ) 


(87 % ) 


(82 %) 


(86 %) 


this work, basis set B 


-0.5276 


-0.3248 


-0.2565 


-0.2230 




(96 % ) 


(94 % ) 


(90 %) 


(97 %) 


exp. (see@) 


-0.555 


-0.345 


-0.285 


-0.229 


LA (see§) 


-0.475 


-0.259 


-0.192 


-0.164 


LDA (seeiS) 


-0.634 


-0.389 


-0.333 




QMC (seel) 


-0.5475 


-0.3587 
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TABLE VI. Lattice constants in Angstrom a) this work, SCF level, basis set A; b) this work, 
core polarization included, basis set A; c) this work, valence correlation included, basis set B - in 
comparison to experimental and other theoretical results. Deviations from experimental values are 
given in parentheses 





C 

V 




p P 


Sn 
Oil 


a) 


3.5590 


5.4993 


5.7516 


6.6001 




(-0.2 To) 


/ I 1 r» 07 \ 

(+1.2 % ) 


(+1.7% J 


( i 1 T07 \ 

(+1.7% J 


b ) 




5.4662 


5.6653 


£} A r* A f\ 

6.4549 






/ i n a 07 \ 
(+0.6 /o ) 


(+0.2/0 J 


( n so/ \ 
(-0.5/0 j 


c) 


3.5833 


5.4256 


5.6413 


6.4443 




(+0.5 %) 


(-0.1 % ) 


(-0.3% ) 


(-0.7% ) 


exp. (seel) 


3.5657 


5.4317 


5.6575 


6.4892 


HF (see0) 


3.58 


5.49 


5.81 




LA (see@) 


3.601 


5.488 


5.760 


6.538 


LDA (see§) 


3.53 


5.38 


5.57 




QMC (seel) 


3.543 


5.404 
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TABLE VII. Bulk moduli (in Mbar) a) this work, SCF level, basis set A; b) this work, core 
polarization included, basis set A; c) this work, valence correlation included, basis set B - in 
comparison to experimental and other theoretical results. Deviations from experimental values are 
given in parentheses 





C 


Si 


Ge 


Sn 


a) 


4.010 


l.Uoo 


U.yol 


U.Ooo 
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FIGURES 



FIG. 1. X n H m clusters treated at the SCF level (H-atoms are not drawn). 

FIG. 2. The X-skeletons of the clusters treated at the CEPA-0 level (big numbers designate 
clusters, small numbers the bonds in each cluster). 
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